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SUMMARY 

This paper describes the use of special interface finite 
elements to model frictional restraint in composite interfaces. 
These elements simulate Coulomb friction at the interface, and are 
incorporated into a standard finite element analysis of a 
two-dimensional isolated fiber pullout test. Various interfacial 
characteristics, such as the distribution of stresses at the 
interface, the extent of slip and delamination, load diffusion 
from fiber to matrix, and the amount of fiber extraction or 
depression are studied for different friction coefficients. The 
results are compared to those obtained analytically using a 
singular integral equation approach , and those obtained by 
assuming a constant interface shear strength. The usefulness of 
these elements in micromechanical modeling of fiber-reinforced 
composite materials is highlighted. 


INTRODUCTION 


In composite materials with brittle constituents, such as 
ceramic matrix composites, the fiber and matrix are unbonded, and 
are restrained by frictional forces at the interface. In such 
composites, the interface has been modeled [1,2,3] as having a 
constant interface shear strength r , and for shear stresses 

greater than t , slipping occurs at the interface. This means of 

analyzing load transfer is approximate, and it is not clear for 
which conditions it is appropriate. 

A more realistic analysis of interfacial load transfer, 
described by Coulomb friction, has been recently presented by 
Dollar and Steif [4]. Instead of assuming a constant interface 
shear strength beyond which slip occurs, the maiximum stress that 
the interface can sustain is related to the prevailing normal 
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stress, c r , and the coefficient of friction at the interface, p, 

n 

by the relation 

x = p |<r I ( 1 ) 

o ^ 1 n 1 

Thus the allowable shear stress changes with the level of the 
applied load, and with the depth of the fiber. 

According to this friction law, at any instant in the loading 
history, either sticking, slipping or opening occurs at a generic 
point along the interface. Conditions for these three states are 
as follows: 

(a) Stick: <r < 0 (compressive) 

n 

M s 


dg/dt = 0, and 

h = dh/dt = 0 (2) 

(b) Slip : o- < 0 

n 

M = mIoJ , 

sgn( dg/dt) = sgn(r), 

h = dh/dt =0 (3) 

(c) Open (equivalent to interface cracking): 

or = 0, 

n 

t = 0, and 

h > 0 (4) 

In these equations, g and h represent the relative slip and 
separation between the fiber and the matrix respectively, while 
d()/dt is the derivative with respect to a time-like parameter 
that increases monotonically sis the loading proceeds. 

The above analysis was used to study the load transfer and 
slip lengths for the simple problem of a two-dimensional fiber 

pullout (and indentation) test (Fig. 1). An analytical integral 

equation approach was used to obtain the solution for the half 
plane problem with the fiber loaded with an axial stress p 
(tension or compression), and subjected to a uniform lateral 
compressive prestress, <r Q . Their results indicate that the extent 

of slipping as well as the rate of load transfer from the fiber to 
the matrix may be signif icantly different from those predicted by 
using the constant shear stress approximation, particularly for 
high compressive stresses, p. In order to avoid added mathematical 
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complexity, the authors have not incorporated the opening case in 
their analysis, and have also assumed that both the fiber and the 
surrounding matrix are made of the same material. 

This paper describes the use of special interface friction 
elements in a finite element model of the problem described. These 
special elements employ incremental elasto-plastic constitutive 
relationships to simulate realistic frictional restraint at 
composite interfaces. This model is capable of implementing all 
the interface conditions mentioned earlier and can also be 
extended to consider non-homogeneous interfaces and more 
complicated finite geometries. 


DESCRIPTION OF THE FRICTION ELEMENT 


A two-dimensional quadratic displacement finite element for 

contact friction problems was developed by Ballarini, Parulekar 

and Plesha [5] for modeling rough, dilatant rock joints, and has 

proved to be quite useful for handling such elastic-plastic 

friction problems. A brief description of the friction model will 

be given here. Each friction element shown in Fig. 2 has zero 

thickness, and 6 nodes with two degrees of freedom at each node 

corresponding to displacements u and u , tangential and normal to 

t n 

the plane of the surface respectively. The constitutive law we 
adopt in this paper was fully developed by Plesha [6], and is 
analogous to the theory of continuum elasto-plasticity. 

The kinematic variables used in the constitutive law are the 
relative surface displacements in the tangential and normal 
directions, which are defined as 

g,= (u B - u 4 )- t 

g = (u - u ) • n , (5) 

n BA 

where u and u are the displacement vectors of adjacent points on 

A B 

the interface associated with the two phases A and B respectively, 
and t and n are unit vectors in the tangent and normal directions 
of the interface respectively. The tangential and normal stresses 
the interface supports at any point are denoted by <r and cr 

t n 

respectively, with the convention that compressive stresses are 
negative. 

A basic assumption in the theory is that the deformation can 
be additively decomposed into 

© . P 
g ■ g, + gj 


i = t, n 


( 6 ) 


4 


where superscripts e and p denote elastic (recoverable) and 
plastic (irrecoverable) parts of the deformation and subscript i 
denotes a vector component in the tangent or normal direction. 
This provides the advantage of leading to a more convenient 
numerical implementation compared to frictional idealizations 
adopted in this paper in which a "stick" condition precedes 
frictional sliding. 

The elastic displacements are related to stresses by the 
elastic constitutive law dtr^s E^dg* , while plastic deformations 
obey the following slip rule 


c 0 if F < 0 or dF < 0 (no slip) 




if F = dF = 0 (slip) 


(7) 


where E are the interface stiffnesses (E , E , E , E ) that 

1 J nn nt 1 1 t n 

have units of force per unit volume, 

gjis the sum of the elastic (pre-slip) and plastic 

(post-slip) displacements at any load level, 

F is a slip function given by F = |<r I + n<r , 

1 t 1 n 

G is a slip potential given by G = |o* t |, and 

X is a non-negative scalar which gives the magnitude of slip. 


For elasto-plastic analysis, the element employs an 
incremental constitutive law relating increments of stress to 
increments of total relative displacements: 


E dg 

{ a<r pq q 

g j~ SF^ 8G~ 

d<r pq d<r 
p q 


dG \ 
dcrj 


( 8 ) 


Since the above explicit relationship is incremental, it is valid 
for arbitrary load and deformation histories that can involve 
loading and subsequent unloading (changes in the direction of 
frictional sliding, for example). 

It has been shown [6] that E = E =0, hence, for no-slip, 

nt t n 

Eq. (8) reduces to 


d <r = E dg , and 
t 1 1 s t 



5 


do- = E dg , or in the matrix form, 

n nn n 



’ d <r 


t 


do- 


ii 



(9) 


while for the slip case, Eq. (8) becomes 
do- = -pE sgn(o- )dg , and 

t nn t n 

do- = E dg , or, 

n nn n 


do- 

t 

_ 

do- 


u n J 


When 

the 


0 -pE sgn(<r ) 

nn t 



’dg t ' 


■ dg n- 


( 10 ) 


fiber and matrix separate under normal tensile 
stress, the interface carries no further stress, Eq. (4). Hence 
for the opening condition, the stiffness of the interface element 
is set to zero. 


Mathematically, E and E are penalty numbers, and in order 
nn 1 1 

to satisfy Eqs. (2)-(4), their values must be manipulated so as to 

obtain the type of rigid-perfect ly plastic response shown in Fig. 

3. Eq. (4) allows no normal separation at the interface until the 

normal stress, <r , turns tensile. Hence E should be as high as 
n nn 

possible for <r < O (compressive), but at the same time should not 
n 

create numerical problems in the solution procedure. Similarly, 
Eq. (2) requires that the elastic part of the tangential 
displacement (g in the pre-slip condition) be invariant with 
respect to the load steps. A parametric study was therefore 
conducted to determine the appropriate stiffness coefficients 
which abide by Eqs. (2) to (4). The following sections present the 
finite element analysis and discussion of the results. 


THE ANALYSIS 


The finite-element model of the half plane fiber pullout (or 
indentation) setup is shown in Fig. 4 (not to scale). Symmetry was 
used to discretize only one half of the domain into 192, 8-node 
isoparametric rectangular elements, while the special 
zero-thickness friction elements were employed at the interface. 
The fiber elements are bound to those of the surrounding matrix 
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through the enforcement of displacement compatibility with the 
interface elements. 

A transient finite element analysis employing an explicit 
time integration scheme developed in Ballarini et. al. [5] is 
implemented in this problem. The load is applied on the fiber in 
small increments, and convergence is observed for step sizes 
greater than 50. The average tangential and normal stresses are 
calculated in each interface element for each load step, and Eqs. 
(2) to (4) are used for determining the stiffness of the interface 
elements at the next load step. 


RESULTS AND DISCUSSION 


In order to check the results of the finite element model, 
various interfacial characteristics, such as the distribution of 
interfacial stresses, the extent of slip and opening, load 
diffusion in the fiber, and slip at the surface are studied for 
various constitutive properties and are compared to those obtained 
by Dollar and Steif [4]. 

Parametric Studies for E and E . 

nn it 

The Young’s modulus, E, for the fiber and matrix is assumed 
as 400 GPa. Different values of E and E were used to study the 

nn 1 1 

behavior of the interface elements with respect to conformity with 

Eqs. (2) to (4). It was observed that in order to satisfy Eqs. (2) 

to (4), the E required is as high as 1.0E15 N/m 3 , while E has 
nn tt 

to be higher than 1.0E14 N/m 3 . 


Interfacial Stresses Under Perfectly Bonded Conditions . 

Perfect bonding between the fiber and matrix is obtained by 
removing the friction elements from the interface. The interfacial 
stress distributions under such conditions was studied in order to 
check the validity of the finite element model, since exact 
solutions for various stresses at the perfectly bonded interface 
are available in the literature. 

The distribution of shear and normal tractions, x and c r 

n 

respectively, due to axial tensile loading of the perfectly bonded 
half-plane with respect to the depth of the fiber is shown in Fig. 
5 for o" 0 =0. These are in excellent agreement with analytical 

solutions for x and <r [4] given by 
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T(d) = ~ 4 P , <r (d) - f f I - tan"‘(d/2) 1 (11) 

ir(d 2 +4) » "l 2 d 2 + 4 J 

where d=depth Into the fiber/fiber radius 
= -y/a 


Variation of Slip Lengths, l g and delamination lengths, 1 Q . 

When the fiber is loaded in tension, the shear stress r is 
maximum at the loading surface (y=0). Slipping takes place when 
the interfacial shear stress t is greater than the allowable shear 
stress r = p I c I. The variation of the extent of slip, 1 , with 

o 1 n 1 s 

load for different friction coefficients, p, is plotted in Fig. 6. 
The extent of slip Increases with load and diminishing friction 
coefficient, since t q is correspondingly smaller. These curves are 

abruptly ended at points where the interface normal stress turns 
tensile, and hence separation occurs at the interface. Fig. 7 
shows the extent of the delamination lengths, 1 for various p. 

The depth of interfacial opening increases with the tensile load 
and diminishing p. 

When the fiber is loaded in compression (indentation), slip 
first initiates at a depth from the surface ranging from one to 
three fiber radii (Fig. 6), and spreads both in depth and towards 
the surface with increasing load. For similar loads, l g for 

compression is less than that for applied tension, particularly 
for high p, where slips initiates at much higher applied loads. No 
opening was observed to occur under compressive loads. The results 
compare very well with those presented in Dollar et. al. [4] (they 
are not duplicated on these plots for the sake of clarity). 


Load Transfer from Fiber to Matrix . 

Stresses are transferred from the fiber to the matrix (or 
vice versa) through the interfacial shear stresses. The variation 
of the axial tensile (and compressive) stresses carried by the 
fiber, <r { , with respect to the fiber depth is shown in Fig. 8. For 

purposes of comparison with the results in Dollar et. al. [4]. the 
applied load was fixed at |p/p<r o |= ± 10, and again, the agreement 

was found to be very good. It is observed that the load diffusion 
is faster for compression than for tension. This is due the fact 
that for compressive loadings, the interface normal stress is 
always compressive, arnd is much higher than that for applied 
tension. Hence the interface is able to sustain much higher shear 
stresses in compression than for tension where extensive slipping 
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and tensile separation occurs thereby restricting the shear 
stresses at the interface . 

The rate of load transfer for a given applied load (p/o- o =±5) 

as a function of p is shown in Fig. 9 (p >0, tensile) and Fig. 10 
(p <0, compressive). As expected, the load diffusion is much 
faster for higher coefficients of friction, and in compression, it 
approaches that for the perfectly bonded interface for 
sufficiently high p (Fig. 10). The load diffusion in the fiber 
under the constant shear stress approximation is shown in the form 
of broken lines in Figs. 9 and 10. The development length, L , 

d 

required in this case to diffuse the fiber load completely is 
given by [4] 

L /a = |p|/r o , (12) 

where t = p<r is the constant shear stress at the interface, 
o o 

The constant interface shear strength approximation neglects 
differences between tensile and compressive loads on the fiber. It 
is observed that for compressive loadings, the diffusion rate is 
faster than that for the constant shear stress approximation, 
while it is slower for applied tension. Unlike the constant shear 
stress approximation, which gives a discrete cutoff load 
dissipation length, L , it is seen that for a frictional 

d 

restraint, a small residual stress is retained in the fiber and 
decreases exponentially with depth. 


Slip at the surface . 

Tensile (or compressive) loads applied on the fiber causes an 
extraction (or indentation) given by the slip between the fiber 
and matrix at the surface. The slip at the surface A , normalized 

S 

by ap(tc + 1)/G (where tc = 3-4v for plane strain, G is the shear 
modulus of the fiber material) is plotted as a function of the 
normalized load p/p<r o in Fig. ll. The results are similar to those 

presented in Dollar et. al. (4], though our results indicate that 
for compression, slip at the surface is somewhat greater than that 
in [4]. For the constant interface shear strength approximation A 

S 

is given by [4] 

A s G/{ap(K+l)> = p/(16po* o ), (13) 

and is plotted in the form of broken lines in Fig. 11. 
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SUMMARY AND CONCLUSIONS 


Preliminary results obtained from the finite element model of 
the pullout test show good agreement with analytical solutions 
predicted by Dollar et. al. [4]. Further investigation is 
necessary to obtain appropriate values for E and E fct for 

inhomogeneous interfaces, and to study stress distributions for a 
system of proximate fibers embedded in the matrix in contrast to 
the case of an isolated fiber. It is seen that interfacial stress 
characteristics for a frictional ly restrained interface is 
different from that obtained from the simple constant shear stress 
approximation, at least for the case of an isolated fiber. It is 
not presently known however whether such frictional modeling would 
affect the overall strength characteristics of composites 
significantly, and whether such a degree of sophistication is 
justified. The interface friction elements may be easily 
incorporated into existing micromechanical analytic models for 
brittle matrix composites (such as the local-global model 
developed by Ballarinl et, al. [7]), which presently assume a 
constant shear strength at the fiber-matrix interface. They could 
also be used to study problems involving complicated and finite 
geometries, and to analyze size effects on interface behavior. 
This would make model assumptions very realistic, and render them 
amenable for the analysis of a wide range of composites. Such 
constitutive local-global modeling using the interface friction 
elements is presently under investigation and will be reported in 
a future communication. 
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Fig. 1. Setup for the fiber pullout test. 
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Fig. 3. The rigid-perfect ly plastic response required of the 
interface element. 
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Fig. 4. Finite element mesh for the fiber pullout test. 

(not to scale) 
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Fig. 5. Distribution of normal (<r ) and tangential (t) interfacial 

a 

stresses for a perfectly bonded interface 
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Fig. 6. Extent of slip as a function of the normalized load p /<r Q . 




1 /a 


17 




18 



Fig. 8. Load diffusion from fiber to matrix for tensile and 
compressive loadings (p/(w Q = ±10) 
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fig- 11. Slip at the surface for tensile and compressive 
as a function of the renormalized load p/po^. 
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